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Abstract 

In this paper, we consider F(R) = R + f(R) theory instead of Einstein gravity with conformal 
anomaly and look for its analytical solutions. Depending on the free parameters, one may obtain 
both uncharged and charged solutions for some classes of F(R) models. Calculation of Kretschmann 
scalar shows that there is a singularity located at r = 0, which the geometry of uncharged (charged) 
solution is corresponding to the Schwarzschild (Reissner-Nordstrom) singularity. Further, we discuss 
the viability of our models in details. We show that these models can be stable depending on their 
parameters and in different epoches of the universe. 



I. INTRODUCTION 

In 1970, Buchdahl considered the first modification of the Einstein Lagrangian density involving an arbitrary 
function of Ricci scalar, F(R) jl|. After that and motivated by inflationary scenarios, Starobinsky proposed an action 
which is to weigh up the effects of R 2 corrections to Einstein gravity |2j. In addition, the recent advent of new 
observational precision tests, such as rotation curves of spiral galaxies Q and solar system tests [H have changed 
the modern view of cosmos based on General Relativity. Amongst the modification of Einstein gravity, the so-called 
F(R) gravity is completely special (see [|| and references therein and for a review see @). Interpretation of the 
hierarchy and singularity problems 0-Q i early-time inflation [lfjj , gravitational waves detection and also the four 
cosmological phases [Hj], have been investigated in F(R) models. 

On the other side, the conformal anomaly (CA) [ljj E3] plays a crucial role in string theory [l5| and gravitational 
field theory (lfj| . Also, it has been shown that there is a deep relation between Hawking radiation and trace anomaly 
[lj} ■ The relation of CA with some cosmological problems such as inflation and cosmological constant problem have 
been considered in Refs. Applications of CA in black hole and particle physics have been studied in fl9l - [2i| 

and [22| . respectively (for a suitable treatment of CA in the black hole solutions see [13] )• In addition, we should note 



that the first Schwarzschild-dS solutions in F(R) gravity has been discussed in [23 1. 

In the semiclassical framework of general relativity, one may use quantum field theory (QFT) to describe the matter 
while the gravitational field is considered as a classical field. When classical gravity is regarded as a background for a 
conformally coupled QFT, the trace of the classical expectation value of matter field energy momentum tensor does 
not vanish. In other word, considering the quantum expectation value of the energy momentum tensor during the 
process of renormalization, leads to additional terms in the Einstein-Hilbert action. 

It has been shown that, when F"(R) ^ 0, metric F(R) gravity is completely equivalent to a special case of 
scalar-tensor theory [24| . In order to find a scalar-tensor representation of F(R) theories, one can use conformal 
transformations [24|. Considering suitable conformal transformations, one may find that, F(R) = R + f(R) behaves 
geometrically as if we have F(R) = R+ (non-minimal scalar field). Remarkably, f(R) plays the role of "scalar field" 
so that the geometry F(R) — R+f(R) becomes equivalent to the Einstein-dilaton geometry in a spherically symmetry 
spacetimc. It is obvious that the CA may come from the quantum effects of matter field. Regarding the scalar-tensor 
representation of F(R) = R + f(R) theories, one may consider f(R) (or equivalently "scalar field") as a matter field. 

It is notable that there are in fact two distinct types of trace anomalies with different geometric and physical 
antecedents. The first anomaly called type A anomaly which is related to the Euler characteristic of 2n-dimensional 
spacetime and the second named B anomaly with a connection to the Weyl tensor. 

In this paper, we consider a class of black hole in F(R) gravity with trace anomaly and look for analytical solutions. 
In order to study the basics solutions of general F(R) gravity theories with constant curvature solutions R = Rq with 
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CA, let us start from the action 

S = Sg + S m , (1) 

in which S g is the 4-dimensional action 

S 9 = T^ fd 4 xV^d[R + f(R)}, (2) 



9 WttG j 

where 8irG = M~ 2 , M p corresponds with the Planck mass, g is determinant of the metric g M „(^, v = 0, 1, 2, 3), R is the 
scalar curvature and R+ f(R) is the function defining the theory under consideration. As the simplest example, the 
Einstein gravity with cosmological constant A is given by f(R) = —2 A. We know that one loop quantum corrections 
lead to a trace anomaly of the energy-momentum (EM) tensor of conformal field theories. In general, as we can show 
that, the trace anomaly has the form [13] 

< >= 9lu/ [l3C aM& C a ^ 5 - a(R 2 - AR aP R a0 + R a0lS R a ^ 5 )] , (3) 

where a, /3 are two positive constants depending on the degrees of freedom of quantum fields and their values are not 
important for our discussions. The first term is a polynomial of Weyl tensor (B anomaly), while the second (Gauss- 
Bonnet) term is Euler characteristic of the 4-dimensional spacetime (A anomaly). From the above considerations, the 
equations of motion (EOM) in the metric formalism are just 

R„ v (1 + fa) - ~ [R + f(R)] g»u + (V M V„ - g^V^fR 

-BnGg^ [(3C a0lS C a ^ s - a(R 2 - 4R aP R a/3 + R aMS R a ^ 5 )] = 0, (4) 

where /r = ^.y . The problem of finding the general static spherically symmetric (SSS) solution in arbitrary f(R) 
theories without imposing the constant curvature condition is in principle too complicated. The required condition 
to get vacuum constant curvature solutions R = i?o(from now i?o will denote a constant curvature value) in vacuum 
implies 

Ro(l + f Ro )=2[R +f(R )] (5) 

For this kind of solutions, an effective cosmological constant may be defined as A e ^ = This type of black- hole 
discussed in [25[. It has been proved that the only static and spherically symmetric vacuum solution (SSVS) with 
constant curvature of any F(R) gravity is just the Hawking- Page black hole in AdS space [26(. But for the theory 
with trace anomaly with T = T£ this condition is 

R (1 + f R ) - 2 [R + f(R)} - 8ttGT = (6) 

One can consider ([6]) as a differential equation (DE) for the f(R) function and with a T — T(r) = T(R). The general 
solution for DE ^ is given by 

f(R) = -R + Cl R 2 + 8nGR 2 f ^j^r-dR' (7) 



II. ASSUMPTION OF EM TENSOR WITH TOPOLOGICAL METRIC 



We consider the external metric for the gravitational field produced by a non rotating object in /(i?)-Weyl (CA) 
gravity. The general 4-dimensional topological metric can be written as [27j 

ds 2 = \{r)dt 2 - li-^dr 2 - r 2 (d6 2 + n 2 d6 2 ) . (8) 

where denotes sin6>i, 1 and sinh6>i for k = 1 (spherical horizon), k = (flat horizon), and k = — 1 (hyperbolic 
horizon), respectively. We follow the methodology of Cai et.al We making an assumption that to have a non 
trivial black hole solution with metric © at a black hole horizon, say r#, the EM tensor must satisfy the auxiliary 
condition, < T t f(r) > r =r H =< 2T(r) > r=rH , which holds not on horizon even in the whole spacetime. We take the 
EM tensor for SSS configurations as 



T£ = diag(p(r), -p(r), -p±{r), -p±{r)), 



(9) 
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where its trace is 



T = p(r) ~p{r) - 2p x - 



(10) 



Applying the auxiliary condition to metric ([8]). one can easily show that A(r) = fi(r). Our main goal is searching 
for an exact solution for (j4|) with definite F(R) action given by ([7]) and determining the metric functions of metric 
([8|). The trace T from Eq. ([T0|) is not constant since the components of the EM tensor (|9]) are functions of the 
radial coordinate r. Thus in general, the f(R) functions are very different depending on the form of the T from (|10|) . 
Locally, let us to take T as a constant which we can find that Eq. (J?} reduces to 



f{R) = -R+ cii? 2 - AttGT. 



(11) 



This action is a second order corrected action of the Einstcin-Hilbert action with a cosmological constant A = 2nGT. 
Thus it's solution is a second order corrected (anti) de Sitter spacetime for a typical (negative) positive value of T. 
Since the Lagrangian of Einstein gravity (i?-term) remove in this model (fTTj) , the corresponding solutions are not 
important for us. In next section we will convert the Eq. ^ to a more convent form. 

III. NEAR HORIZON AND APPROXIMATE SOLUTIONS FOR CONSTANT CURVATURE 



In this section we focuss on the special case f(R) — — 2A, with constant curvature i.e. R = Rq. Firstly it is 
adequate to rewrite the EM ([3]) as the following form 



where we used the following identity [2 



•- - a ) R 2 - 2 08 - 2a) R a pR af) + (fi - a)R aMS R a ^ s ) 



(12) 



ot/3y6 



l 



-R z — 2R a pR a ^ + Rap^sR 



(13) 



As we pointed former, a, /? are two constants depending on the degrees of freedom of quantum fields and their values 
are not important for our discussion. We choice a = j3 and also define a new parameter a' = 2a. Thus, we write the 
next EOM 



1 



8irGa' 



iWi + fR)-[^[R + f(R)] - ^p# 2 



+ (V M V„ - g^VpV^fR = 87TGa'g^(R a pR a0 ). 
Constant curvature solutions for ([T^| with f(R) = — 2A, can be written as 

Rn» =u"g^{j + R a pR a0 ), 



(14) 



(15) 



where 7 = 



and a" = 8irGa' . Applying the metric ([5]) with A = fi to Eq. (|15p . one leads to the next 



nonlinear differential equations (NDEs) 



A"(r) 



-2rX'(r) ± W2r 4 



-7 + 



8a" 



4[rA'(r) + A(r) - k] 



1 

2a 77 



-2rA'(r) ± W2r 



-7 + 



-rA'(r) +k - A(r) 



-4[r\'{r) + X(r)-k} 2 



(16) 



(17) 



where prime denotes the derivative with respect to the radial coordinate r. We must solve Eqs. (fl~6| and (fl7|) 
simultaneously. This NDEs have no (simple) exact solution but we can solve it approximately or numerically with a 
suitable boundary conditions. According to the Hawking-Bekenstein formula for temperature of the black hole, if the 
metric ^ posses a black hole solution with a horizon located at r — h, we can deduce X(h) =0 and 



T 



47T 



(18) 
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Examining (|16|) and (|17l) at horizon we can write 



X"(h) = 



-8nhT± , 2ft 4 



-7 + 



8n 



irl 



-A{AirhT-kf 



1 

2^" 



-8nhT± , 2h 4 



-7 + 



-AirhT + k 



t"h? 



- A(AnhT - fc) 2 



(19) 



(20) 



Here, we desire to write near horizon solution of A(r). To do this, we equate Eq. (|16p with Eq. (|17[) and obtain 



3/i 2 [SirhT - 2 Ah 2 + R h 2 - 2k] 



12 [AnhT - fc) - R 2 h 4 



(21) 



In other word, Eq. (|16p is equal to Eq. (|17p for specific value of a" , (I2"TT) . Thus, the near horizon solution of A(r) of 
Eqs. ([H]) and flU]) with Eq. (J2U) is obtained as 



A(r) ~ 4vrT(r - /i) + 



(r - hf 
2b? 



-8ichT± J2h 4 



-7 + 



-4nhT + k 

a c h 2 



4(4irhT - kf 



(22) 



On the other side, for small values of a", we can expand Eqs. (|16l) and (fl~7)) and solve the resulting DEs with the 
following solution 



1M r 2 

X( r )=k + -(2A-R ). 

r b 

Identifying this metric function with a Schwarzschild-(a)ds metric gives 

Ro = 4A, 

and thus, for small values of a", the metric ([5]) reads as 



ds 2 = 



, 2M r 2 , k 

k + — 2A - R 

r 6 



dt 2 



dr 2 



[fc-2M + f (2A-.Ro)] 



r 2 dfl 2 



(23) 



(24) 



(25) 



IV. EXACT SOLUTIONS 
A. f(R) = -2A model: 
1. Case I: a ± 0, P = 

At the first step and following the procedure in [2l[, we consider f(R) = — 2A with type A anomaly (Gauss-Bonnet 
term). It is easy to show that for A(r) = fi(r), the field equation ^ reduces to 

4- [n + ry = o, (26) 

ar 

where 

E = 2a"X'(r) [A(r) - A] + (27) 

^ = i ^A'(r) i = M ' (28) 
Using Eq. ([26)) with Eqs. (|27|) and (l28l) . the metric function can be obtained as 

X(r) = k-Ar 2 , (29) 
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where 

^ 3 ± v^W (30) 

We should note that for a" = 3/(16A), Eq. (|30[) reduces to 2A/3 and the presented solution is close to asymptotic 
(a)dS solution. 



2. Case III: arbitrary a = ft 

In what follows, we obtain an exact solution of the field equation (|16p with A and B anomalies, simultaneously 
with f(R) — — 2A, which is Einstein-A gravity with CA. Straightforward calculations show that, in this model, one 
can obtain topological Schwarzschild solution as [23[ 

9M 

A(r) =fi(r) =k-Ar 2 -— , (31) 

r 

where A is given in Eq. (f30]) and j3" — 16ttG(3 — a". As we mentioned before, one can obtain approximate 
Schwarzschild- (a) dS solution, provided the parameters of the solution are chosen suitably. 



B. f(R) = dR 2 - 2A model: 

Considering a linear function of T(R') = cR' + b, we find that Eq. (0 leads to the Einstein-Hilbert action with a 
R 2 correction with a cosmological constant A = 2bnG and c = . Straightforward calculations show that Eq. (f3"Tj) 
is a solution of field equation Q for a" — ft" and obtained A (|3Df. 

We note that asymptotic flat Schwarzschild solution may be found by setting fc = 1 and 6 = 0, which exactly yields 
the vanishing cosmological constant. In other word, R 2 correction does not effect on the asymptotic behavior of the 
solutions. 

In addition to the Schwarzschild solution, one can show that topological charged solution is another exact solution 
of this model. In other word, we should note that 

A(r) =Ai(r) = k-Ar 2 + ^, (32) 
is a solution of FE (|4]) , provided the parameters of the solution are chosen as follow 

P" = (33) 

Cl - 241' ^ 

where A is the same as Eq. (|30"|) . 

It is notable that the presented uncharged black holes is the same as solutions which obtained in Ref. [3(| for a = 0. 
Although for g tt — g~J~ in Ref. |30| , the metric function has a charged term pj- in analogy with our charged solution, 
but these solutions arc completely different. They have different geometry and Ricci scalar R with various asymptotic 
behavior. 

Here, we want to discuss about the stability of such solution as an example of F(R) gravity solutions with CA. The 
stability of the solutions in F(R) gravity has been discussed in literatures [3l|. In fact, the stability of the de Sitter 
solution may be obtained by imposing the one-loop effective action to be real. In order to study of stability, some 
various techniques have been employed which all these methods are in agreement with the following conditions which 
state the existence and the stability of the de Sitter solution: 

/ : 2^) = R F'(R ), 

11 ^! >0 

• F'(Ro) >U ' 

HI : F ' {Ro) - > 1. (35) 
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Equations I and II, state the existence of a solution with positive constant curvature, while equation 777 ensures the 
stability of such a solution. We should note that one can obtain such conditions by a classical perturbation method 

& 

Here we use these conditions to our de Sitter solutions. In this model one has 

F(R) = R-2A + Cl R 2 , 
F'(R) = l + 2 Cl R, 
F"(R) = 2ci. 

Substitute de Sitter solution in which i?o = 4 A (A > 0), we obtain 

F(R Q ) = 2A+16ciA 2 , 
F'(Rq) = l + 8ciA, 
F"(Ro) = 2c x . 

which confirm that the first and second conditions of Eq. (|35[) satisfied and for third condition we get 

1+ 8^A >1 ' 

As a result, de Sitter solution of this model is stable for c\ > 0. 



C. f(R) = Cl R 2 + Kln(R) - 2A model: 

The cosmological evolution of the f(R) models based on logarithmic correction has been studied and it has been 
claimed that such theories have a well-defined Newtonian limit |33| . Here, we desire to consider R 2 with logarithmic 
corrections, simultaneously. To do this, one can consider T(R') = cR' + ailn(R') + a2 in Eq. ([7]), which leads to 
f(R) = c x R 2 +K\n(R) - 2A, when 

' (36) 



8irG 

* - -«r (37) 

K + 4A 

Considering the mentioned f(R), one can show that Eq. (|3"Tj) is a solution of field equation for 

13" = a", (39) 
4(A + 12AV--3A) 

21n(12A)-l ' 1 ' 

It is notable that for k = 1, one may found asymptotic dS solution for specific value of conformal parameter a" = 
^[21n(4A)-l]. 

Now, we search about charged solution. Considering the mentioned model with Eq. (j40|) . it is easy to show that 
Eq. (EH) with P" = 5a"/6 satisfy the FE flU), provided 

_ -1 A + 12a" A 2 + 6A [ln(12A) - 1] 
Cl_ 72 A 2 (21n(12A)-l) ' ( ^ 



D. f(R) = Cl R 2 + c 2 R n - 2A model with n > 2: 

Although most of f(R) gravity models are well-behavior in the weak gravity regime (34[, But for the strong gravity 
regime, some of these models have a serious drawback such as singularity problem. In Ref. (35| . Kobayashi and Maeda 
has been resolved the singularity problem arising in the strong gravity by considering a higher curvature correction 
term proportional to R n . In order to obtain the mentioned f(R) model (R 2 and R n corrections to Einstein- Hilbert 
action), one may consider T{R') = ^R' + ( " 8 ~ R ' n + 2§U in Ec l- ©■ 
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To obtain Schwarzschild solution Eq. (131) in this model, it is sufficient to consider p = a" with 

-A{k+12A 2 a" -ZA) 
° 2 (n-2)(12A)" ' 1 ; 

which confirm that depend on the sign of A, one may achieve asymptotic adS/dS solutions provided a" = — 3c2(n — 
2)(4A)"" 2 . 

Considering /?" = 5a"/6 with Eq.f |4"2"|) . it is straightforward to obtain charged solution Eq. (1521 for 

1 nA + Una" A 2 - 6A(n - 1) 
Cl= 72 AHn-2) (43) 

We should note that the former solution is asymptotically dS (adS) for A > (A < 0) if we set the free parameters 
of this model in the following manner 

1 ,A AS 

Cl = -8A + 6(^)< (44) 
a" 

° 2 ~ _ 3(n-2)(4A)™- 2 ' ( } 

V. PROPERTIES OF THE EXACT SOLUTIONS 

Now, let us look for the singularity of the solutions. Considering Schwarzschild- (a)dS like (|3 1 1) with metric <JSj> , one 
can show that, the Kretschmann scalar is given by 

R^R^ = 2AA 2 + , (46) 

therefore there is a singularity located at r = (with nonzero a — p), as the Schwarzschild solution in general 
relativity. 

In addition, it is easy to show that in the case of charged solution (|3"2"|) , the Kretschmann scalar is 

R^ KV R^ V = 2AA 2 + -|-, (47) 

which confirm that the geometry of singularity (r — > 0) is the same as Reissner-Nordstrom black hole. It is notable 
that presented charged black hole solution (l3"2"j) , has two horizons for hyperbolic horizon (fc = —1) and for k = 0, 1, 
we encounter with naked singularity (see Fig. [T]for more details). 

It is notable that the presented solutions are nonsingular for vanishing mass (M) and charge (q). 

VI. VIABILITY CONDITIONS OF F(R) THEORIES 

The first stud y on cosmological viable f(R) models have been presented in Ref. (36{. The fundamental conditions 
and restrictions |37| that are usually imposed to action 

S « = Teb / d ^F{R) (48) 
theories to provide consistent both gravitational and cosmological evolutions are: 

1. F"(R) > for R ^> F"(R). This is the stability requirement for a high curvature classical regime [38[ and 
that of the existence of a matter dominated era in cosmological evolution. A simple physical interpretation 
can be given to this condition: if an effective gravitational constant G e // = G/(l + F'(R)) is defined, then 
the sign of its variation with respect to R, dG e ff/dR, is uniquely determined by the sign of F"(R), so in case 
F"{R) < 0, G e ff would grow as R does, because R generates more and more curvature itself. This mechanism 
would destabilize the theory, as it wouldn't have a fundamental state because any small curvature would grow 
to infinite. Instead, if F"(R) > 0, a counter reaction mechanism operates to compensate this R growth and 
stabilize the system. 



8 




FIG. 1: A(r) (Eq. ([32}) versus r for A = -0.1, g = 1, and k = 1 (bold line), k = (solid line) and fc = -1 (dashed line). 



2. 1 + F'(R) > 0. This conditions ensures that the effective gravitational constant is positive, as it can be checked 
from the previous definition of G e ff- It can also be seen from a quantum point of view as the condition that 
avoids the graviton from becoming a ghost (39j . 

3. F'(R) < 0. Keeping in mind the strong restrictions of Big Bang nucleosynthesis and cosmic microwave back- 
ground, this condition ensures GR behavior to be recovered at early times, that is, F(R)/R — > and F'(R) — > 
as R — > oo. Conditions 1 and 2 together demand F(R) to be a monotone increasing function between the values 
-KF'(i?)<0. 

4. F'(R) must be small in recent epochs. This condition is mandatory in order to satisfy imposed restrictions by 
local (solar and galactic) gravity tests. As the analysis done in indicates, the value of \F'(R)\ must not be 
bigger than 10 -6 (although there is still some controversy about this). This is not a needed requirement if the 
only goal is to obtain a model that explains cosmic acceleration. 

Now we examine the viability conditions (1-4) for a general F(R) = R + f(R) model described in (7). Explicitly, 
we have 

F'(R) = 2c x R + WnGR [" ^^-dR' + 8nG^^ (49) 
J R' A R 

F"{R) = 2ci + 8ttG^^- + 16nG /"* ^^-dR' (50) 
R z J /£ M 

Validity of the conditions (1-4) depends on the form of the T(R). We will summarize the different models of 
F(R)=R+f(R) discussed in this paper as the follows. Note that since we considered the F(R) = R + f(R) mod- 
els, thus in each case we firstly written the explicit for of the F(R). 



A. T(R)=T=c.t.e 

In this case, one can obtain 



F'(R) = 2 Cl R (51) 
F"(R) = 2ci (52) 



with the following properties: 
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1. F"(R) — 2c\R > for R ^> 2c\. Thus this model is stable or a high curvature classical regime. The effective 
gravitational constant is G e // = \ + 2 Ci r - ^ n tne eaj ly universe epoch, when R » R = 2c\, this effective G 

reduces to G e // ~ ,its variation with respect to R, d G d ^ f < 0. Thus G e // would decreases as R grows in the 
early times, this mechanism would stabilize the theory, as it would have a fundamental state . 

2. Obviously for R >• 2c\,l + F'(R) > 0. It means that the G e ff is positive. Thus this case of theory is free from 
ghosts. 

3. F(R)/R — > oo and F'(R) — > oo as R — > oo. Thus only for ci = 0, this model recovers GR behavior at early 
times. 

4. In early times, F'(R) «. Thus our model passes the local solar system tests. Further, from |2cii?| < 10~ 6 gives 
the bound on |ci| < ^-,R ~ 0(A). 



B. f(R) = -2A 



This case describes the Einstein-Hilbert action with a cosmological constant with 

F'(R) = 1 (53) 
F"(R) = (54) 

1. F"(R) = 1 > for ReR. Thus this model is stable . 

2. Obviously 1 + F'(R) > 0. It means that the G e ff — G is positive. Thus this case of theory is free from ghosts. 
Further, it coincides with GR at the action level. 



C. f(R) = Cl R 2 - 2A 



Now, we have F(R) = R + c x R 2 - 2A and thus we've 

F'(R) = 1 + 2 Cl R (55) 
F"(R) = 2ci (56) 

1. F"(R) = 2c\ > for c\ > 0. This is the stability requirement for a high curvature classical regime of this 
model and that of the existence of a matter dominated era in cosmological evolution. This model induces an 
effective gravitational constant G e fj = G/2(l + c\R) is defined, then the sign of its variation with respect to R, 

~lm L = ~ 2{i+cIr)' 1 ' ^ s uniquely determined by the sign of ci, so in case c\ < 0, G e // would grow as R does, 
because R generates more and more curvature itself. The problem changed when the sign of the c\ reverses. This 
mechanism would destabilize the theory, as it wouldn't have a fundamental state because any small curvature 
would grow to infinite. Instead, if ci > 0, a counter reaction mechanism operates to compensate this R growth 
and stabilize the system. 

2. 1 + F'(R) = 2(1 + ciR) > 0. This conditions ensures that the effective gravitational constant is positive and 
from a quantum point of view as the condition for avoiding from a ghost graviton. 

3. 2c\R < for c\ < O.This condition ensures GR behavior to be recovered at early times. 

4. F'(R) = 1 + 2c\R must be small in recent epochs. This condition is mandatory in order to satisfy imposed 
restrictions by local (solar and galactic) gravity tests. The value of |1 + 2c\R\ must not be bigger than 10~ 6 . 

D. f(R) = Cl R 2 + K ln(i?) - 2A 

In this case the Lagrangian of the model is F(R) = R + c\R? + Kh\(R) — 2A, we have 

F'{R) = l + 2 Cl R+^- (57) 
R 

F"(R)=2 Cl -^ (58) 
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1. F"(R) — 2ci — J^z > for R > J This is the stability requirement for a high curvature classical regime of 

this model and that of the existence of a matter dominated era in cosmological evolution. This model induces an 
effective gravitational constant G e // = G/(2 + 2c\R + 4|) is defined, then the sign of its variation with respect 

to R, d = — (2*4.20 1 R+i-)' 2 ' * s UIU( l ue ly determined by the sign of c%, K, so in case c\ < 0, K > 0, G e // would 
grow as R does, because R generates more and more curvature itself. The problem changed when the sign of the 
ci,K reverses. This mechanism would destabilize the theory, as it wouldn't have a fundamental state because 
any small curvature would grow to infinite. 

2. 1 + F'(R) = 2 + 2cxR + f > 0. This conditions ensures that the effective gravitational constant is positive and 
from a quantum point of view as the condition for avoiding from a ghost graviton. The root of the equation 
1 + F'(R) = 2 + 2c x R + § = locates at R± = ~ 1±v ^~ 8ciK -. This is the critical curvature. If ci > then the 
model is free from ghosts for R > R + , R < R_ . But if Cj < then the ghosts are present for this range of the 
curvature. 

3. 1 + 2c\R + j£ < for c\ < ,R > R +} R < i?_ . This condition ensures GR behavior to be recovered at early 
times. 

4. F'(R) = 1 + 2c\R + must be small in recent epochs only for small negative values of c\. This condition 
is mandatory in order to satisfy imposed restrictions by local (solar and galactic) gravity tests. The value of 
|1 + 2c x R + -|| must not be bigger than 10~ 6 . 



E. Model with f(R) = Cl R 2 + c 2 R n - 2A with n > 2 

In this case, we have F(R) = R + c\R 2 + C2R 71 — 2A. Explicitly, we have 

F'{R) = 1 + 2 Cl R + nczR 71 - 1 (59) 
F"(R) = 2ci + n(n — l)c 2 -R"~ 2 (60) 

Validity of the conditions (1-4) depends on the signs of the constants Ci, c 2 for n > 2. If we take C\ > 0, c 2 > 0, then 
obviously for all values of the R,F'(R) > 0,F"(R) > 0. The conditions (1,2) are satisfied automatically. But for 
c\ < it is possible to have F'(R) < . Thus the model satisfy the local solar system tests for some values of the 
C\,C2 ■ For example take n=3. Then the condition F'(R) < can be satisfied by such values of the curvature in the 
interval i?_ < R < R + for c 2 > 0,cf > 3c 2 and R > R + or R < J?_ for c 2 < where R± = —c\ ± y/cf — 3c 2 . 



F. On the validity of CA coefficients 



Anomaly is a classical symmetry breaking by quantum corrections. It seems that by such a definition, this effect 
has significance role only in early universe. But now, we show that in a universe which is in accelerating expansion 
phase, filled by dark energy, the phenomena can be described by CA terms. The main idea belongs to [41| . We 
summarized here the main results (for details refer to the 4l| and references there). As we know that, the CA term, 
if it is considered as an IR contribution to the stress-energy, helps to understanding of behavior of the EM tensor of 
the matter fields at low energy of any effective theory. As we know that, the trace term < > is of order of the 

< >^ H 4 where H is the usual Hubble constant. Cut-off from QCD tells us that this value must be of order 

< Tjf >~ 10" 12 (e^) 4 42]. This is the observational vacuum energy density. Difference between these values is the 
first cosmological constant problem and has no satisfactory solution till now. By a systematic analysis of the CA 
terms corrections to the vacuum energy density it has been shown that under certain conditions the CA terms can 
provide an appropriate scale observed in nature as predicted by QCD estimates. The first step is writing the total 
action for the low energy theory which contains the usual Einstein-Hilbert action, the conformally invariant Weyl 
portion, the higher order curvature invariants and finally the anomalous contribution. By writing the total EM tensor 
we observe that the resulting expression is correctly proportional to the equations of motion for some auxiliary fields. 
It is easy to show that the anomalous contribution to the stress tensor and the Weyl one would be of the same order of 
magnitude. If the spacetime been conformally flat, then the Weyl invariant term vanishes, we have only the anomalous 
contribution as the only unique contribution to the energy tensor. Consider the spacetime as the static de Sitter. By 
calculating the finite contribution to the energy and pressure density, it is shown that the result does not depend on 
the CA coefficients. Thus we remove the divergences of the CA stress tensor. Indeed, the contribution of the CA 
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terms reduces to the local, second order in the curvature, geometrical terms, which are of order H 4 , and therefore 
irrelevant at scales lower than the Planck one. 



VII. CONCLUSION 

It has been shown that F(R) = R + f(R) gravity has (a)dS Schwarzschild solution. In addition, charged solutions 
of some pure F(R) models have been investigated in Rcf. [5]. In the present paper, we considered CA as a source of 
F(R) gravity and investigated its consequences. 

At first, we used temperature conception of black holes to find a near horizon solution and then we found an 
approximate solution which one might interpret it as a Schwarzschild- (a) dS solution. 

Next step is devoted to find exact solutions. Calculations showed that depending on the values of free parameters, 
one may find Schwarzschild and charged solutions. Also, we found that the geometry of uncharged (charged) singularity 
(r — > 0) is corresponding to the Schwarzschild (Reissner-Nordstrdm) black hole. 

Finally, we investigated the fundamental conditions to obtain viable F(R) gravity. In other word, we found some 
restrictions on free parameters to guarantee stability as well as ghosts free conditions. In addition, we checked the 
behavior of the F(R) theory at early times and also satisfaction of local (solar and galactic) gravity tests. 

It is worthwhile to investigate the thermodynamics properties and also generalize our solutions to higher dimensional 
spacetime, and these problems are left for the future. 
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